AD-A270  768 

IllllllilP 


HEAT  TRANSFER  PERFORMANCE  OF  A  ROOF-SPRAY 
COOLING  SYSTEM  EMPLOYING  THE 
TRANSFER  FUNCTION  METHOD 


— » 


By 

JOSEPH  A.  CLEMENTS 


This  document  has  been  approved 
for  public  release  and  sols;  its 


A  THESIS  PRESENTED  TO  THE  GRADUATE  SCHOOL 
OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
MASTER  OF  SCIENCE 

UNIVERSITY  OF  FLORIDA 


*^-<f 


93-24488 

niiliiiii 


1993 


This  thesis  is  dedicated  to  my  wife,  Maria,  whose  patience  and 
endurance  made  this  thesis  possible. 


ACKNOWLEDGMENTS 


I  wish  to  express  my  sincere  thanks  to  Dr.  S.  A.  Sherif,  Dr.  H.  A.  Ingley  III,  and 
Dr.  B.  L.  Capehart  for  consenting  to  serve  as  members  of  my  supervisory  committee.  The 
guidance,  patience,  and  encouragement  of  Dr.  Sherif,  my  committee  chairman,  have  been 
greatly  appreciated. 

My  graduate  education  has  been  funded  by  the  United  States  Navy's  Fully  Funded 
Graduate  Education  Program.  I  am  most  grateful  to  Captain  J.  P.  Phelan  for  his  generous 
and  skilled  assistance  in  dealing  with  the  Navy  that  was  the  key  to  my  being  able  to 
maintain  my  educational  funding  and  remain  in  Gainesville  full-time  to  complete  my 


degree. 


TABLE  OF  CONTENTS 

pagg 

ACKNOWLEDGMENTS  iii 

LIST  OF  TABLES  .  vi 

LIST  OF  FIGURES  .  vii 

KEY  TO  SYMBOLS  . viii 

ABSTRACT  .  x 

CHAPTERS 

1  INTRODUCTION  .  1 

2  MATHEMATICAL  MODEL  . 4 

Introduction  . 4 

Heat  Transfer  Mechanisms  . 4 

Roof  Exterior-Air  Interface  . 5 

Roof  Inside  Surface-Room  Air  Interface  . 8 

Thermal  Capacitance  Effects  of  Roof  . 9 

Combined  Mathematical  Model  .  11 

Experimental  Roof  Parameters  .  12 

Weather  Model  and  Data  .  12 

3  RESULTS  .  17 

Introduction  .  17 

Roof  Heat  Transfer  and  Temperatures  at  B*-y an,  Texas  .  17 

Roof  Heat  Transfer  at  Pittsburgh  .  24 

Summary  .  26 


IV 


4 


CONCLUSIONS  AND  RECOMMENDATIONS 


29 


Conclusions  .  29 

Recommendations  .  31 

APPENDICES 

A  COMPUTER  CODE  LISTING  33 

B  UNCERTAINTY  ANALYSIS  39 

REFERENCES  .  42 

BIOGRAPHICAL  SKETCH  .  44 


V 


LIST  OF  TABLES 


TABLE  PAGE 

2.1  Thermal  Characteristics  of  Roof  at  Bryan,  Texas  .  14 

2.2  Thermal  Characteristics  of  50  mm  Thick  Concrete  Roof  at 

Pittsburgh  .  15 

2.3  Thermal  Characteristics  of  50  mm  Thick  Pine  Plank  Roof 

at  Pittsburgh  .  16 

3.1  Experimental  and  Predicted  Heat  Flux  and  Temperature 
Comparison  for  Roof-Spray  Cooled  Roof  at  Bryan, 

Texas  .  21 

3.2  Experimental  and  Predicted  Weather  for  Bryan,  Texas, 

and  Pittsburgh  .  22 

3.3  Experimental  and  Predicted  Heat  Flux  and  Temperature 

Comparison  for  Dry  Roof  at  Bryan,  Texas  .  24 

3.4  Experimental  and  Predicted  Heat  Flux  Comparison  for 

Roof-Spray  Cooled  Concrete  Roof  at  Pittsburgh  .  28 

3.5  Experimental  and  Predicted  Heat  Flux  Comparison  for 

Roof-Spray  Cooled  Pine  Plank  Roof  at  Pittsburgh  .  28 

B.l  Uncertainty  Analysis  by  Sequential  Perturbation  for 

Concrete  Roof  at  Pittsburgh  .  40 


LIST  OF  FIGURES 


FIGURE  PAGE 

2.1  Energy  Balances  at  Roof  Surfaces  . 6 

3.1  Relation  Between  Time  and  Heat  Flux  Trough  Outside 

and  Inside  Surfaces  of  Roof-Spray  Cooled  Roof  at 

Bryan,  Texas,  in  July  .  20 

3.2  Relation  Between  Time  and  Temperature  of  Outside 
Surface  of  Roof-Spray  Cooled  Roof  at  Bryan,  Texas,  in 

July  .  20 

3.3  Relation  Between  Time  and  Heat  Flux  Trough  Outside 

and  Inside  Surfaces  of  Dry  Roof  at  Bryan,  Texas  in  July  ...  23 

3  .4  Relation  Between  Time  and  Temperature  of  Outside 

Surface  of  Dry  Roof  at  Bryan,  Texas,  in  July  .  23 

3.5  Relation  Between  Time  and  Heat  Flux  Through  Outside 

and  Inside  Surfaces  of  Concrete  Roof  at  Pittsburgh  .  27 

3.6  Relation  Between  Time  and  Heat  Flux  Through  Outside 

and  Inside  Surfaces  of  Pine  Plank  Roof  at  Pittsburgh  .  27 


KEY  TO  SYMBOLS 


A  Area  of  roof, 

A,  Constants  in  water  saturation  pressure/temperature  correlation 

equation 

b.  Conduction  transfer  function  coefficients,  W/(m^  ®C) 

c„  Conduction  transfer  function  coefficients,  W/fm^  ®C) 

d.  Conduction  transfer  function  coefficients 

Evaporative  heat  transfer  coefficient,  5.678  W/(m^  K) 
hj  Internal  convective  heat  transfer  coefficient,  W/(m^  K) 

h„  External  convective  heat  transfer  coefficient,  W/(in^  K) 

n  Summation  index 

P  Total  barometric  pressure  of  moist  air.  Pa 

P„  Partial  pressure  of  water  vapor  in  moist  air.  Pa 

P„  Saturation  pressure  of  water  at  the  roof  surface.  Pa 

P„  Saturation  pressure  of  water  vapor,  psia  or  Pa 

flco.<i(fooo  Conductive  heat  flux  through  roof,  W/m^ 

Qco«v(iiuide)  Convective  heat  flux  at  inside  surface  of  the  roof,  W/m^ 

flco.v(oBt.ide)  Convective  heat  flux  at  outside  surface  of  the  roof,  W/m^ 
qe,tf  Heat  gain  through  roof  at  calculation  hour  $,  W/m^ 

^«v.p  Evaporative  heat  flux,  W/m^ 


vui 


^r.d(i..ide)  Radiative  heat  flux  at  inside  surface  of  the  roof,  W/m^ 
^r»d(o.uide)  Radiative  heat  flux  at  outside  surface  of  the  roof,  W/m’ 


Qtour  Solar  heat  flux  at  outside  surface  of  the  roof,  W/m^ 

R,  Thermal  resistance  of  roof  (m^  K/W) 

T  Temperature,  °R  or  K 

T,  Ambient  outside  temperature,  “C  or  K 

T^  Sol-air  temperature,  “C  or  K 

T^  i  Temperature  of  inside  surface  of  the  roof,  “C  or  K 

T,  „  Temperature  of  outside  surface  of  the  roof,  ®C  or  K 

T„,  Temperature  of  inside  conditioned  space,  K 

V  Wind  speed,  k/hr 

W  Uncertainty 

Wx,  Uncertainty  of  any  variable  Xj 

a  Absorptance  of  the  surface  of  the  roof  for  solar  radiation 

6  Time  interval,  h 

5R  Difference  between  the  long-wave  radiation  incident  on  the  surface 

from  the  sky  and  surroundings  and  the  radiation  emitted  by  a 

blackbody  at  ambient  outdoor  air  temperature,  W/m^ 

8  Emissivity 

a  Stefan-Boltzman  constant,  5.6697  x  10  *  W/(m^  K"*) 

$  Time,  h 

0}  Humidity  ratio  of  the  moist  air,  kg„„ykgj^ 


IX 


Abstract  of  Thesis  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Master  of  Science 

HEAT  TRANSFER  PERFORMANCE  OF  A  ROOF-SPRAY 
COOLING  SYSTEM  EMPLOYING  THE 
TRANSFER  FUNCTION  METHOD 

By 

Joseph  A.  Clements 
December,  1993 


Chairperson:  Dr.  S.  A.  Sherif 

Major  Department.  Mechanical  Engineering 

Roof-spray  cooling  systems  have  been  developed  and  implemented  to 
reduce  the  heat  gain  through  roofs  so  that  conventional  cooling  systems  can  be 
reduced  in  size  or  eliminated.  Currently,  roof-spray  systems  are  achieving 
greater  effectiveness  due  to  the  availability  of  direct  digital  controls.  The 
objective  of  this  thesis  was  to  develop  a  mathematical  model  of  the  heat 
transfer  through  a  roof-spray  cooled  roof  that  predicted  heat  transfer  based 
on  existing  weather  data  and  roof  heat  transfer  characteristics  as  described  by 
the  Transfer  Function  Method.  The  predicted  results  of  this  model  were 
compared  to  the  results  of  existing  experimental  data  from  previously 
conducted  roof-spray  cooling  experiments. 


X 


The  mathematical  model  is  based  on  energy  balances  at  the  exterior  and 
interior  surfaces  of  the  roof  construction  that  include  evaporative,  convective, 
radiative,  and  conductive  heat  transfer  mechanisms.  The  Transfer  Function 
Method  is  used  to  relate  the  energy  balances  at  the  two  surfaces  that  differ  in 
amplitude  and  phase  due  to  the  thermal  resistance  and  thermal  capacitance 
characteristics  of  the  roof. 

The  combined  model  yielded  moderately  good  predictions  of  heat 
transfer  through  the  roof  exferimental  results  for  the  roof-spray  cooled 
condition.  The  calculation  method  shows  promise  as  a  relatively  simple  means 
of  predicting  heat  gains  based  on  calculation  procedures  that  are  similar  to 
those  frequently  used  by  many  engineers. 
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CHAPTER  1 
INTRODUCTION 


Roof-spray  cooling  systems  continue  to  increase  in  popularity  for  a 
number  of  reasons.  Evaporative  roof-spray  cooling  which  has  been  widely 
implemented  since  the  1930s  to  provide  cost  effective  cooling  for  industrial 
applications  is  now  seen  as  a  method  to  provide  cooling  which  does  not  add  to 
the  depletion  of  fossil  fuels  or  contribute  to  possible  global  warming  problems 
caused  by  combustion  and  the  use  of  CFC  based  refrigerants  in  traditional 
vapor  compression  refrigeration  systems.  A  major  advance  that  has  improved 
the  effectiveness  of  roof-spray  cooling  is  the  availability  of  direct  digital 
controls,  which  allow  water  to  be  applied  to  roofs  to  achieve  maximum 
cooling  without  the  application  of  excessive  water. 

Many  researchers  have  studied  various  forms  of  rooftop  evaporative 
cooling  including  solar  ponds,  trickling  water,  and  roof-sprays.  In  1940 
Houghton  et  al.  [1]  studied  the  effects  of  roof  ponds  and  roof-sprays  on 
temperature  and  heat  transfer  in  various  roof  constructions.  In  the  1960s, 
Yellott  [2]  reported  on  the  effectiveness  of  a  roof-spray  cooling  system  that 
consisted  of  little  more  than  a  rooftop  grid  of  sprinkler  heads  and  supply  pipes 
with  a  water  supply  controlled  by  a  solenoid  valve  and  timer.  More  recently. 
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research  has  been  conducted  to  numerically  model  roof-spray  cooling  systems. 
This  includes  the  work  of  Tiwari  et  al.  [3],  Carrasco  et  al.  [4],  Somasundaram 
and  Carrasco  [S],  and  Kondepudi  [6]. 

The  heat  transfer  through  roofs  with  various  construction 
configurations  has  also  been  studied  extensively.  Spitler  and  McQuiston  [7] 
discussed  the  latest  modeling  techniques  for  the  prediction  of  cooling  loads 
and  times  as  they  occur  in  the  zones  of  a  building  as  a  result  of  the  daily  loads 
of  ambient  temperature  and  incident  solar  radiation. 

However,  a  model  that  incorporates  the  current  modeling  of  roof 
systems  and  roof-spray  cooling  is  lacking.  This  thesis  presents  a  numerical 
model  for  the  combined  heat  transfer  response  to  roof-spray  cooling  and  roof 
construction  in  response  to  the  variation  of  ambient  temperature,  humidity  and 
solar  flux  over  a  typical  day  and  incorporates  the  effects  of  thermal 
capacitance  of  a  roof  as  modeled  by  the  Transfer  Function  Method  (TFM).  The 
TFM  uses  conduction  transfer  functions  to  predict  the  heat  flux  at  the  inside 
of  the  roof  as  a  function  of  previous  values  of  outside  and  inside  temperatures 
and  heat  fluxes.  The  magnitude  and  direction  of  the  inside  heat  flux  may  be 
smaller  or  larger  than  the  heat  flux  at  the  outside  surface  depending  on 
whether  the  thermal  mass  of  the  roof  is  increasing  or  decreasing  in  stored 
energy.  This  thesis  will  present  a  numerical  solution  to  the  heat  transfer 
problem  and  compare  the  results  to  existing  experimental  data  from  previously 
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conducted  roof-spray  cooling  experiments.  The  heat  transfer  mechanisms 
acting  on  the  outside  surface  in  response  to  ambient  temperature  and  incident 
solar  radiation  are  evaporation,  convection,  radiation,  and  conduction.  The 
transfer  mechanisms  acting  on  the  inside  surface  are  convection,  radiation,  and 
conduction  in  response  to  the  energy  transfer  at  the  outside  surface  and  a 
fixed  room  ambient  temperature.  The  energy  balances  at  the  two  surfaces  are 
related,  but  the  amplitude  and  phase  of  the  inside  energy  transfer  will  depend 
on  the  thermal  resistance  and  thermal  capacitance  characteristics  of  the  roof. 

The  solution  presented  should  enhance  the  understanding  of  the 
combined  effects  of  evaporative  cooling  and  thermal  capacitance  of  building 
envelope  construction.  As  a  design  tool  the  benefit  of  this  modeling  will  be 
that  building  designers  would  more  easily  choose  the  best  combination  of 
roofing,  spray-cooling,  and  conventional  cooling  systems.  This  modeling  also 
has  an  added  benefit  in  that  it  can  be  used  as  a  basis  for  algorithms  that 
determine  the  most  effective  and  efficient  application  of  roof-spray  cooling  to 
a  roof 
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CHAPTER  2 

MATHEMATICAL  MODEL 

Introduction 

The  various  elements  of  the  heat  transfer  mechanisms  are  presented  in 
this  chapter.  First,  the  mathematical  basis  (derived  in  terms  of  the  energy 
balances)  is  presented  in  the  form  of  equations.  The  different  heat  transfer 
mechanisms  and  assumptions  concerning  these  mechanisms  are  also  discussed. 
Second,  the  thermal  capacitance  characteristics  and  effects  of  the  roof  system 
are  presented.  The  different  types  of  roof  construction,  the  heat  transfer 
properties  of  the  materials,  and  the  assumptions  concerning  these  properties 
are  also  discussed.  Finally  the  interaction  of  the  heat  transfer  mechanisms  and 
the  construction  are  discussed  and  a  mathematical  basis  in  the  form  of 
equations  is  presented. 

Heat  Transfer  Mechanisms 

Kondepudi  [6]  and  Carrasco  et  al.  [4]  presented  a  thorough 
development  of  a  mathematical  model  of  the  heat  transfer  mechanisms  of  a 
spray  cooled  roof  with  some  assumptions  being  made.  These  assumptions 
include; 
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(1)  The  problem  is  quasi-steady  in  response  to  varying  ambient 
temperatures  and  solar  flux. 

(2)  An  infinitesimally  thick  film  of  water  is  maintained  on  the  roof 
surface. 

(3)  The  saturated  water  vapor  pressures  respond  in  a  linear  fashion 
to  temperature. 

(4)  Wind  effects  are  included  in  the  convective  heat  transfer  model. 

(5)  Thermal  capacitance  effects  of  the  roof  and  water  film  are 
neglected. 

All  of  these  assumptions  with  the  important  exception  of  the  thermal 
capacitance  effects  of  the  roof  will  be  carried  through  this  paper.  The  initial 
formulation  of  the  energy  balances  will  ignore  the  thermal  capacitance  effects 
of  the  roof 

The  energy  balances  for  the  two  interfaces  of  the  process  are  presented 
in  Figure  2.1.  The  energy  balances  at  each  interface  are  described  below  and 
provide  a  problem  that  must  be  solved  simultaneously. 

Roof  Exterior-Air  Interface 

As  indicated  in  Figure  2.1  the  incident  solar  energy  is  dissipated  by  a 
combination  of  evaporative,  radiative,  convective,  and  conductive  heat 
transfer. 


Qtoltr  flev»p'^Qfi«l(o»l»ide)'^flco«v(o«liid«)'^^co«d<toof) 
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Figure  2. 1 :  Energy  balances  at  roof  surfaces 


The  evaporative  heat  transfer  is  a  result  of  Pick's  Law  of  division  of  the  water 
vapor  at  the  wetted  surface  diffusing  into  the  ambient  air,  as  presented  by 
Tiwari  et  al.  [3]. 


q.v.p  0  013h„.,(P„-P„)  (2.2) 

The  partial  pressure  of  water  at  the  roof  surface,  P„,  is  predicted  from 
Mathur's  correlation  [8]  for  the  water  saturation  pressure,  P,„  as  a  function  of 
temperature. 


P„  =Ao  +  A,T  +  AJ'  +  A3T'  +  AJ^ 


where  Aq  =  670.2012  psia 
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A,  =-5.325521  psia/“R 
A,  =  1.59464  xlO-2  psia/“R' 

Aj  =-0.2134061  X10-^  psia/“R" 

A^  =  0.1077853  xlO"’  psia/^R^ 

P„  =  saturation  pressure  of  water  vapor,  psia 

T  =  absolute  temperature,  "R 

The  partial  pressure  of  water  in  the  ambient  air  is  predicted  by  ASHRAE  [9] 
as; 


P„  =  (Pa})/(0.62198  +  «) 


where  P  =  total  barometric  pressure  of  moist  air.  Pa 
CO  =  humidity  ratio  of  the  moist  air, 

9r»d(o«ii<«)~®®^('rr,c"Tj)  (2.3) 

The  equation  for  th<-  outside  radiative  heat  flux  can  be  refined  based  on 
Alford  et  al.  [10],  who  proposed  that  the  ambient  temperature  for  the  radiative 
heat  transfer  should  be  decreased  by  2.8  K.  This  decrease  is  based  on  the 
assumption  that  all  radiation  reaching  the  earth  originates  in  the  first  1500 
meters  of  elevation  and  accounts  for  a  normal  atmospheric  temperature 
gradient  of  about  7  K  per  1000  meters  elevation.  Thus, 

q,.d(o.««.re‘^(T?,o-<T.-2.8)^)  (2.4) 

Another  method  of  correcting  for  the  difference  in  net  longwave 
radiation  exchange  between  the  roof  and  the  sky  is  presented  by  Bliss  [11]. 
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Under  clear  sky  conditions  the  sky-to-surface  radiation  difference  for  a 
horizontal  black  surface,  5R,  is  typically  in  the  range  of  63  to  95  W/m^  This 

radiation  difference  is  a  function  of  elevation,  cloudiness,  air  temperature  near 
the  ground,  and  the  dew  point,  but  will  be  assumed  to  be  an  average  of  63 
W/m^  for  clear  sky  conditions  at  the  experiment  sites  of  Bryan,  Texas  and 
Pittsburgh,  Pennsylvania.  Bliss'  model  will  be  used  in  the  predicted  heat  flux 
calculations  in  this  report. 

The  convective  and  conductive  heat  fluxes  at  the  outside  roof  surface 
are  given  by  the  following  two  equations: 


4<aav(oatti4«) 


(2.5) 


9coad(rooO  (^r.o“Tf,i)/^r 


(2.6) 


Roof  Inside  Surface-Room  Air  Interface 

As  indicated  in  Figure  2.1,  the  heat  transfer  mechanisms  acting  at  the 
inside  surface  are  conduction,  radiation,  and  convection. 

4coad(toof)  ^i«l(ia*id«)~^^c<iav(ia<id«)  (2.7) 


^MdCiBiide)  ®^(Tr,i”Tset) 


(2.8) 


9coav(iiiiide)  ^i('^t,i”'^set) 


(2.9) 
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On  the  basis  of  weather  data  from  a  database  that  includes  solar 
radiation,  ambient  temperature,  and  humidity,  the  above  equations  can  be 
solved  iteratively  to  obtain  a  predicted  roof  surface  temperature  and  a 
predicted  conductive  heat  flux  into  the  roof.  This  is  based  on  Kondepudi's 
assumption  [6]  that  there  are  no  thermal  capacitance  effects  of  the  roof  or 
water  film 

Thermal  Capacitance  Effects  of  Roof 

The  above  discussion  provides  a  simple  methodology  for  predicting  the 
outside  roof  temperatures  and  heat  fluxes  if  the  thermal  capacitance  effects  of 
the  roof  are  neglected.  In  actual  applications  the  thermal  mass  effects  of  the 
building  envelope  can  play  an  important  role  in  decreasing  and/or  shifting  the 
phase  of  the  cooling  load  on  a  building.  Several  load  calculation 
methodologies  have  been  developed  to  account  for  this  thermal  mass  effect. 
Spitler  and  McQuiston  [12]  presented  a  review  of  the  most  popular  cooling 
load  calculation  methodologies  that  are  being  used  for  manual  and  computer 
calculations.  Because  of  the  proliferation  of  computers  in  calculating  cooling 
loads.  The  American  Society  of  Heating,  Refrigerating,  and  Air-Conditioning 
Engineers  (ASHRAE)  [13]  is  currently  concentrating  its  research  on  improving 
and  expanding  the  range  of  circumstances  to  which  the  TFM  can  be  applied. 

The  TFM,  as  described  by  Mitalas  [14],  was  developed  in  the  late 
1960's  as  a  computer-oriented  load  calculation  method  that  applies  weighting 
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factors  to  the  loads  experienced  by  the  exterior  surfaces  of  a  building  to 
account  for  the  thermal  storage  effects  of  the  building  envelope.  These 
weighting  factors  are  referred  to  as  conduction  transfer  function  (CTF) 

coefficients  and  relate  the  difference  in  the  sol-air  temperature  and  the  inside 
space  temperature  to  predict  the  difference  in  magnitude  and  phase  shift 
between  the  load  on  the  outside  surface  of  the  structure  and  the  load  that  is 
converted  to  a  cooling  load.  Sol-air  temperature  [13]  is  an  equivalent  outside 
air  temperature  that,  in  the  absence  of  all  radiation  changes,  accounts  for  the 
heat  transferred  to  the  exterior  surface  due  to  the  combination  of  incident 
solar  radiation  and  radiative  and  convective  heat  transfer  between  the  exterior 
surface  and  the  ambient  air  based  on 

T  =T +aq„,„ /h,-e6IUio  (2.10) 

which  is  derived  from  the  following  two  equations  [13]; 

q=«q...„ +h,(T.-T,>e5R  (2.11) 

q=h„(T.-T,,)  (2.12) 

Based  on  the  sol-air  temperature  described  above  and  a  constant  inside  room 
air  temperature,  the  TFM  predicts  the  heat  gain  through  a  roof  by  the 
following  equation. 


The  CTF  coefficients  take  into  account  the  heat  transfer  coefficients  at  the 


exterior  and  interior  surfaces  of  the  roof  and  the  roof  construction.  CTF 
coefficients  can  be  calculated  using  a  computer  program  presented  by 
McQuiston  and  Spitler  [7]  or  use  approximate  values  of  CTF  coefficients  that 
were  developed  based  on  investigations  by  Harris  and  McQuiston  [IS]  into  the 
thermal  characteristics  of  walls  and  roofs  and  their  effect  on  the  heat  transfer 
to  conditioned  spaces.  These  approximate  values  are  available  in  tabular  form 
in  various  texts  (see  [7]  and  [13]  ,  for  instance)  and  are  included  in  many 
computerized  load  calculation  programs  such  as  the  software  provided  in 

Carrasco  et  al.  [4]  and  the  Load  Design  Program  [16]. 

Combined  Mathematical  Model 

The  models  for  evaporative  roof-spray  cooling  by  Kondepudi  [6]  and 
the  Transfer  Function  Method  (TFM)  described  above  were  combined  to 
predict  the  effect  of  evaporative  roof-spray  cooling  and  the  thermal 
capacitance  effect  of  roof  construction.  Due  to  the  iterative  nature  of  the 
solution  of  the  roof-spray  cooling  problem  and  the  application  of  the  TFM,  the 
use  of  a  computer  to  solve  the  combined  problem  is  virtually  imperative.  The 
problem  is  solved  in  a  two-step  process  that  first  ignores  the  thermal 
capacitance  effect  of  the  roof  to  predict  a  heat  flux  into  the  roof  and  then 
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applies  the  TFM  to  predict  the  cooling  load  on  the  interior  space  that  results 
from  the  calculated  exterior  heat  flux  This  solution  predicts  the  amplitude 
decay  in  the  load  as  a  result  of  the  evaporative  cooling  and  the  roof 
construction  and  the  phase  shift  in  the  cooling  load  as  a  result  of  the  thermal 
behavior  of  the  roof  in  response  to  imposed  loads. 

Experimental  Roof  Parameters 

The  thermal  characteristics  of  the  roof  that  Somasundaram  et  al  [5] 
used  in  their  experiments  are  provided  in  Table  2.1,  which  includes  conduction 
transfer  coefficients  obtained  from  Reference  [16],  The  thermal  properties  of 
each  layer  were  taken  from  measurements  made  by  Somasundaram  et  ai.  [5] 
and  from  tabular  data  from  Reference  [13].  A  study  of  the  resulting  amplitude 
decay  and  phase  shift  in  the  cooling  load  is  conducted  in  Chapter  3. 

Weather  Model  and  Data 

The  weather  data  that  was  used  for  the  prediction  of  heat  transfer  in  the 
Bryan,  Texas  roof  was  taken  from  the  Typical  Meteorological  Year  (TMY) 
readings  for  Waco,  Texas  averaged  over  the  period  of  the  experimental 
roof-spray  cooling.  The  data  of  Waco,  Texas  were  used  because  it  is  located 
only  75  miles  from  Bryan,  Texas,  which  has  similar  summer  design 
temperatures,  and  it  is  only  70  meters  higher  in  elevation  [13].  The  ASHRAE 
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solar  heat  flux  calculation  method  [13]  which  is  based  on  the  work  by 
Threlkeld  and  Jordan  [17]  could  also  be  used  TMY  weather  tape  was 
selected  for  its  solar  radiation  data  that  makes  it  a  good  choice  for  solar 
design  problems. 

Actual  experiment  site  weather  measurements  including  ambient  wet- 
and  dry-bulb  temperatures,  solar  radiation,  and  wind  speed  were  reported  by 
Houghten  et  al.  [1]  for  the  experiments  in  Pittsburgh.  Ambient  pressure  data 
were  extracted  from  Trane  Typical  Cooling  Day  weather  files  for  Pittsburgh 
[16].  Typical  Cooling  Day  pressure  varied  only  340  Pa,  so  a  constant  pressure 
could  have  been  assumed  with  little  effect  on  predicted  results. 
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Table  2  .1:  Thermal  Characteristics  of  Roof  at  Bryan,  Texas 


Physical  Description 


Description 

Thickness  Conductivity 
{mm}  {W/m^F} 

Density 

{kg/m^} 

Specific  Heat 
{kJ/(kgK)} 

Resistance 
{(m^  KVW} 

Outside  air  film 

0 

000 

0 

00 

044 

Rigid  insulation 

25 

.043 

32 

.84 

587 

Felt  and  asphalt 

10 

.192 

1121 

1  68 

.050 

Gypsum  roof  deck 

150 

156 

800 

1  09 

978 

Mineral  insulation 

150 

069 

10 

.84 

2.201 

Insulation  air  film 

0 

.000 

0 

00 

.120 

Decrement  factor  =120,  Overall  weight  =  135  kg/m^ ,  Overall  U-value  =  250  W/(m^  K) 


Time  delay  =  1 1  Hours, 

Heat  capacity  =  337  kJ/(m^  kg  K),  C,  coefficient  =  183E-02 

n  B.  coefficients 

D,  coefficients 

0  .1733146E-08 

lOOOOOOE+01 

1  .1402007E-04 

-.2086841E-H)1 

2  .3365204E-03 

.1422010E-K)1 

3  .9128446E-03 

-.3563488E+00 

4  .5027138E-03 

.2923254E-01 

5  .63786 16E-04 

-.7691044E-03 

6  .1774670E-05 

.5973486E-05 
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Table  2.2.  Thermal  Characteristics  of  50  mm  Thick  Concrete  Roof  at 
Pittsburgh 


Physical  Description 


Description 

Thickness  Conductivity 
{mm}  {W/m^  K} 

Density 

{kg/m^ 

Specific  Heat 
{kJ/(kgK)} 

Resistance 
{(m'  KVW} 

Outside  air  film 

0 

000 

0 

00 

,044 

Slag 

13 

1.436 

880 

1.68 

009 

Felt  and  asphalt 

10 

.192 

1120 

1.68 

050 

Slag 

13 

1.436 

880 

1.68 

009 

Felt  and  Asphalt 

19 

.192 

1120 

1.68 

.050 

Concrete 

50 

1.731 

2240 

.84 

.029 

Inside  air  film 

0 

.000 

0 

,00 

.120 

Decrement  factor  =  .775,  Overall  weight  =  168  kg/m^ ,  Overall  U-value  =  .346  W/(m^  K) 
Time  delay  =  4  Hours,  Heat  capacity  =  410  kJ/(m^  kg  K),  C,  coefficient  =  653 


n  B,  coefficients  D,  coefficients 

0  .3732012E-01  lOOOOOOE+Ol 

1  .3805050E+01  -.9006944E+00 

2  .2269784E-K)0  .1360358E+00 

3  .9101764E-02  -.2234380E-03 

4  .7930927E-05  .5431017E-08 

5  .8646022E-11  0 

6  0  0 
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Table  2.3:  Thermal  Characteristics  of  50  mm  Thick  Pine  Plank  Roof  at 
Pittsburgh 


Physical  Description 


Description 

Thickness  Conductivity 

Density 

Specific  Heat 

Resistance 

{mm} 

{W/m^  F} 

{kg/m^} 

{kJ/(kgK)} 

{(m^  KyW} 

Outside  air  film 

0 

000 

0 

00 

044 

Slag 

13 

1.436 

880 

1.68 

Felt  and  asphalt 

10 

192 

1120 

1.68 

Slag 

13 

1.436 

880 

1  68 

Felt  and  Asphalt 

19 

192 

1120 

1.68 

Pine  Planking 

50 

.121 

590 

2.51 

419 

Inside  air  film 

0 

.000 

0 

.00 

.120 

Decrement  factor  =  708,  Overall  weight  =  84  kg/m^ 

,  Overall  U-value  =  1  33  W/(m'  K) 

Time  delay  =  5  Hours, 

Heat  capacity 

=  370  kJ/(m' kg  K),  C. 

coefficient  = 

189 

n  B,  coefficients 


D,  coefficients 

o 


0  .3732012E-01 

1  .3805050E+01 

2  .2269784E-K)0 

3  .9101764E-02 

4  .7930927E-05 

5  .8646022E-11 

6  0 


lOOOOOOE+Ol 

-.9006944E+00 

.1360358E-K)0 

-.2234380E-03 

.5431017E-08 

0 

0 


CHAPTER  3 
RESULTS 

Introduction 

This  section  presents  the  numerical  results  of  predicted  cooling  loads 
along  with  the  experimental  results.  First,  the  predicted  results  for  a  structure 
in  Bryan,  Texas  are  compared  to  the  experimental  results  obtained  by 
Somasundaram  and  Carrasco  [S].  Then  predicted  cooling  loads  for  two  roofs 
at  Pittsburgh  are  compared  to  experimental  results  of  Houghten  et  al.  [1], 

Roof  Heat  Transfer  and  Temperatures  at  Br\  ttn.  Texas 

The  comparisons  of  the  experimental  data  and  numerical  predictions  for 

the  roof-spray  cooled  roof  at  Bryan,  Texas  are  plotted  in  Figures  3.1  and  3.2 

and  provided  in  Table  3.1.  The  predicted  average  heat  gain  through  the  roof, 
q,,  of -4.7  ±1.2  W/m^  is  in  very  good  agreement  with  the  experimental  value 
of  -4.2  W/m^.  The  predicted  outside  roof  surface  temperatures  are  in  good 

agreement  with  the  experimental  data,  with  an  average  overprediction  of  less 

than  one  "C.  The  peak  predicted  conductive  heat  flux  through  the  outside  of 

the  roof,  qc„,d(,ooO’  lower  than  the  peak  experimental  data.  The  average 

predicted  qe„.d(,ooo  **  90%  higher  than  the  experimental  results.  Some 
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variations  are  expected  due  to  the  thermal  capacitance  effects  of  the  roof  that 
were  neglected  in  the  prediction  of  the  outer  conductive  heat  flux  and  the  fact 
that  the  outside  flux  meter  was  located  below  the  gypsum  roof  deck  in  the 
experiments. 

The  results  plotted  in  Figures  3.1  and  3.2  also  indicate  two  likely 
problem  areas  with  the  models  used  to  predict  results.  These  are  weather  and 
the  thermal  properties  of  the  roof.  First,  although  the  cumulative  daily  solar 
radiation  obtained  from  the  TMY  weather  data  [17]  is  very  near  the 
experimental  amount.  Table  3.2,  the  peak  experimental  solar  radiation  was 
likely  higher  at  midday  and  tapered  off  more  rapidly  at  times  other  than 
midday.  This  would  account,  in  part,  for  the  di^erences  in  the  shapes  of  the 
predicted  and  experimental  outside  surface  heat  flux  and  temperature  curves. 
The  second  likely  weather  problem  was  an  underprediction  of  the  net 
nocturnal  radiation  exchange  between  the  roof  and  the  atmosphere  that  was 
reflected  in  the  underprediction  of  the  nighttime  cooling  effect.  This  also 
would  have  contributed  to  the  predicted  average  q^iddoco  being  90%  higher 
than  the  experimental  average  value.  Minor  inaccuracies  in  the  thermal 
properties  of  the  roof  model  including  the  convective  heat  transfer  coefficients 
also  would  have  contributed  to  some  of  the  differences  in  the  experimental  and 
predicted  results.  This  can  be  illustrated  by  the  fact  that  the  predicted 
difference  between  the  maximum  and  minimum  values  of  q,  is  0.7  W/m^  while 
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the  experimental  difference  is  1.6  W/m^  which  may  be  a  result  of  inaccuracies 
in  the  mass  and  thermal  conductivity  of  the  roof  materials  used  in  the 
mathematical  model  of  the  roof  provided  in  Table  2.1. 

As  a  check  on  the  validity  of  the  roof  and  exterior  load  models, 
experimental  results  from  dry,  no  roof-spray  cooling,  tests  on  the  same  roof 
during  the  period  of  17-21  August  were  compared  to  predicted  results.  The 
comparison  of  experimental  data  and  numerical  predictions  for  the  dry  roof  are 
plotted  in  Figures  3.3  and  3.4  Due  to  the  21%  difference  in  experimental  and 
TMY  cumulative  solar  radiation  data,  the  TMY  hourly  solar  radiation  data 
were  adjusted  upward  to  match  the  experimental  daily  cumulative  data.  The 
TMY  hourly  ambient  temperatures  averaged  only  1.1  ®C  lower  than  the 
experimental  data  and  were  not  adjusted. 

The  predicted  difference  between  the  maximum  and  minimum  values  of 

q,  is  1.2  W/m^  while  the  experimental  difference  is  1.3  W/m^  which  indicates  a 

good  representation  of  the  thermal  capacitance  effects  of  the  roof  The 
predicted  average  q,  of  2.4  ±1  W/m^  is  in  moderately  good  agreement  with 
the  experimental  value  of  0.8  W/m^.  This  difference  may  be  a  result  of  a 

combination  of  inaccuracies  in  the  weather,  heat  transfer  model,  and  thermal 

characteristics  of  the  roof  that  are  magnified  by  the  higher  heat  fluxes 

associated  with  the  dry  roof  condition. 


ROOF  SURFACE  TEMPERATURE  (DBG  Q 


SOUR  FLUX  ON  ROOF  SURFACE  (W/m») 
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Table  3.1:  Experimental  and  Predicted  Heat  Flux  and  Temperature 
Comparison  for  Roof-Spray  Cooled  Roof  at  Bryan,  Texas 


July  17  -  July  30,  1987 
Outside  surface  Inside  surface 


Experimental 

Predicted 

Experimental 

Predicted 

Average  Roof  Heat  Flux 
(W/m') 

0.40 

0.76 

0.44 

0.40 

±0.2 

Maximum  and  Minimum 

7.1 

3  9 

1.3 

0.8 

Roof  Heat  Flux 
(W/m') 

-3.1 

-2.1 

-0.2 

0.1 

Cumulative  Daily  Roof 
Heat  Flux 
(W  h/m') 

9.5 

18.3 

10.5 

10.0 

Average  Roof  Surface 
Temperature 

rc) 

28.2 

29.1 

26.7 

Note:  1 

Maximum  and  Minimum 

41.4 

41.2 

27.7 

Note:  1 

Roof  Surface 
Temperature 
("C) 

20.3 

18.8 

25.8 

Note:  1 

Note;  1.  Experimental  data  for  inside  surface  temperature  were  used  in  the 
predict  other  table  values. 


Source; 


Experimental  data  from  reference  [5]. 
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Table  3.2  Experimental  and  Predicted  Weather  for  Bryan,  Texas,  and 
Pittsburgh 


Bryan,  Texas 

Pittsburgh 

Wet  Roof 

Dry  Roof 

Wet  Roof 

Test  Date 

Test  Date 

Test  Date 

July  17- 

August  17- 

8  September- 

July  30 

August  21 

9  September 

Experimental  Predicted 

Experimental  Predicted 

Experimental 

Cumulative  Daily 

6398  6497 

6492  5114 

3613 

Solar  Radiation 

Note;  1 

(W  hr/m^  /24  hrs) 

Average  Ambient 

28.54  28.64 

30.24  29.12 

26.3 

Temperature 

CO 

Maximum  and  Minimum  34.79  34.44 

37.49  33.89 

35.0 

Ambient  Temperatures 

23.12  22.72 

24.33  25.00 

18.6 

CO 

Average  Wind  Speed 

Note;  2  14.2 

Note:  2  13.6 

22.3 

(k/hr) 

Maximum  and  Nfinimum  Note;  2  17.8 

Note:  2  16.7 

38.6 

Houriy  Average 

Note:  2  10.2 

Note;  2  10.6 

9.7 

Wind  Speed 
(k/hr) 


Notes:  1 .  Due  to  large  difference  between  predicted  and  experimental  values, 
solar  data  were  adjusted  to  experimental  value  for  predicting  other 
table  values. 

2.  No  experimental  wind  data  available. 

Sources;  Experimental  data  from  Reference  [S]. 

Predicted  data  from  Reference  [17]. 


ROOF  SURFACE  TEMPERATURE  (DBG  Q 


0  2  4  6  8  10  12  14  16  18  20  22  24 

TIME  OF  DAY  (HR) 

Figure  3.4:  Relation  Between  Time  and  Temperature  of  Outside  Surface  of 
Dry  Roof  at  Bryan,  Texas,  in  July 
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Table  3.3.  Experimental  and  Predicted  Heat  Flux  and  Temperature 
Comparison  for  Dry  Roof  at  Bryan,  Texas 


August  17  -  August  21,  1987 


Outside  surface 

Inside  surface 

Experimental  Predicted 

Experimental 

Predicted 

Average  Roof  Heat  Flux 

3.4 

2.4 

0.8 

2.4 

(W/m') 

±10 

Maximum  and  Minimum 

17.6 

9.7 

1.5 

3.0 

Roof  Heat  Flux 
(W/m^) 

-4.1 

-1.4 

0.2 

18 

Cumulative  Daily  Roof 
Heat  Flux 
(W  h/m') 

82.5 

58.5 

19.3 

58.2 

Average  Roof  Surface 

Temperature 

(“C) 

38.8 

37.2 

28.3 

Note! 

Maximum  and  Minimum 

66.7 

63.6 

29.4 

Note;  1 

Roof  Surface 

21.1 

23.1 

26.9 

Note;  1 

Temperature  (®C) 


Note;  1.  Experimental  data  for  inside  surface  temperature  was  used  in  the 
prediction  of  other  table  values. 

Source;  Experimental  data  from  reference  [5]. 


Roof  Heat  Transfer  at  Pittsburgh 


The  numerical  predictions  and  experimental  data  for  the  concrete  roof, 
which  are  given  in  Table  3.4  and  in  Figure  3.S,  show  good  agreement.  No 
experimental  outside  roof  heat  flux  or  temperature  data  are  available. 
However,  the  experimental  and  predicted  inside  heat  flux  for  the  concrete  roof 
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are  in  good  agreement  in  hourly  magnitude,  phase  shift,  and  cumulative  daily 
heat  flux  into  the  room.  The  predicted  average  q,  of  -4  7  ±12  W/m^  is  in  very 
good  agreement  with  the  experimental  value  of  -4  2  W/m^  The  predicted 

difference  between  the  maximum  and  minimum  values  of  q,  is  28  :  W/m^  34% 

higher  than  the  experimental  difference,  21.1  W/m^  The  differences  in  the 

predicted  and  experimental  inside  heat  fluxes  may,  in  part,  be  attributed  to 

minor  inaccuracies  in  the  mathematical  descriptions  of  the  roof  and  the  models 

of  heat  transfer  at  the  roof  surfaces. 

The  numerical  predictions  and  experimental  data  for  the  pine  plank 

roof,  which  are  given  in  Table  3.S  and  in  Figure  3.6,  show  moderately  good 
agreement.  The  average  predicted  q,  of  -4.2  ±1.4  W/m^  is  in  moderately  good 
agreement  with  the  experimental  value  of  -2.6  W/m^  The  difference  between 

the  predicted  and  experimental  values  and  the  approximate  two  hour  difference 

in  the  phases  of  the  predicted  and  experimental  q,  indicate  that  the 

mathematical  description  of  the  pine  plank  roof  is  inaccurate.  The 

mathematical  descriptions  of  the  Pittsburgh  roofs  are  based  on  sketches  by 

Houghten  et  al.  [1]  that  provide  no  thermal  characteristics  of  the  roof 

materials.  There  is  a  large  uncertainty  in  the  thermal  response  modeling  of  the 

pine  plank  roof  because  of  the  unknown  characteristics  of  the  pine  planking. 

It  is  unknown  if  the  two  inch  thickness  of  the  pine  planking  given  by  Houghten 

is  the  actual  or  nominal  thickness.  Also,  the  density,  thermal  conductivity, 

and  specific  heat  of  the  wood  could  vary  by  over  30%  depending  on  the 
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specific  species  and  moisture  content  of  the  planking.  An  example  of  the 
variation  in  the  thermal  properties  of  a  material  from  expected  is  the  thermal 
conductivity  of  the  Bryan,  Texas  mineral  insulation  that  testing  [S]  revealed 
was  60%  higher  than  values  found  in  Reference  [9] 

Summary 

The  computer  code,  based  on  the  mathematical  models,  provides  a 
moderate  to  good  comparison  with  experimental  results  of  heat  flux  through 
the  inside  surface  of  a  roof  and  into  a  room.  The  predicted  average  heat  gain 
through  the  roof  for  the  Bryan,  Texas  insulated  concrete  roof,  the  Pittsburgh 
concrete  roof,  and  Pittsburgh  pine  plank  roof  are  in  good  agreement  with  the 
experimental  averages.  The  prediction  of  heat  flux  and  temperature  at  the 
outside  surface  of  the  roof  provided  good  qualitative  comparison  with 
experimental  results.  However,  large  magnitude  and  small  phase  differences 
from  experimental  data  were  obtained  as  a  result  of  the  thermal  capacitance 
effects  of  the  roof  that  were  neglected  in  the  prediction  of  the  outside  surface 
temperature  and  heat  flux  and  the  placement  of  the  outside  roof  heat  flux 
meter  below  the  gypsum  decking. 
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Figure  3.5;  Relation  Between  Time  and  Heat  Flux  Through  Outside  and 
Inside  Surfaces  of  Concrete  Roof  at  Pittsburgh 


Figure  3.6:  Relation  Between  Time  and  Heat  Flux  Through  Outside  and 
Inside  Surfaces  of  Pine  Plank  Roof  at  Pittsburgh 
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Table  3  4  Experimental  and  Predicted  Heat  Flux  Comparison  for 
Roof-Spray  Cooled  Concrete  Roof  at  Pittsburgh 


Outside  surface 
Predicted 

August  8,  1939 

Inside  surface 
Experimental  Predicted 

Average  Roof  Heat  Flux 
(W/m') 

10 

-4.2 

-4.7 

±1.2 

Maximum  and  Minimum 
Roof  Heat  Flux 
(W/m') 

15.6 

-12.9 

6.3 

-14.8 

9  1 
-19  1 

Cumulative  Daily  Roof 
Heat  Flux 
(W  h/m2) 

24.6 

-99  5 

-113.5 

Source:  Experimental  data  from  reference  [  1  ] 

Table  3.5:  Experimental  and  Predicted  Heat  Flux  Comparison  for 
Roof-Spray  Cooled  Pine  Plank  Roof  at  Pittsburgh 

August  8,  1939 

Outside  surface  Inside  surface 

Predicted  Experimental  Predicted 

Average  Roof  Heat  Flux 
(W/m^) 

0.9 

-2.6 

-4.2 

±1.4 

Maximum  and  Minimum 
Roof  Heat  Flux 
(W/m^) 

16.4 

-11.9 

2.4 

-8.5 

7.1 

-15.9 

Cumulative  Daily  Roof 

Heat  Flux 

(Wh/m^) 

21.2 

-62.3 

101.7 

Source:  Experimental  data  from  reference  [1]. 


CHAPTER  4 

CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

The  objective  of  this  thesis  was  to  combine  existing  models  for 
evaporative  roof-spray  cooling  and  the  transient  heat  flow  through  roofs  to 
provide  a  model  that  will  predict  roof  cooling  loads  based  on  known  weather 
data  and  roof  material  descriptions.  The  best  evaluation  of  the  effectiveness 
of  the  model  is  by  comparison  to  well-documented  experimental  data. 

Two  roof  conditions  were  studied.  One  was  dry  roof  (no  spray 
cooling),  and  the  other  was  damp  roof  (spray  cooled).  Three  different  roof 
constructions  at  two  locations  were  studied.  The  roofs  included  an  insulated 
gypsum  deck,  an  uninsulated  concrete  deck,  and  an  uninsulated  pine  plank 
deck. 

The  mathematical  formulation  gave  a  simple  but  effective  prediction  of 
the  roof  cooling  load  for  an  interior  that  is  maintained  at  a  fairly  constant 
temperature.  The  agreement  of  experimental  data  with  predictions  of  roof 
cooling  loads  as  a  result  of  the  combined  effects  of  roof-spray  cooling  and 
roof  construction  gave  a  high  level  of  confldence  that  the  use  of  the  combined 
roof-spray  cooling  and  Transfer  Function  Method  model  will  yield  good 
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predictions  for  buildings  maintained  at  a  constant  internal  temperature  The 
agreement  between  predicted  and  experimental  average  heat  gain  through  the 
roof,  q,,  values  was  very  good  making  the  predicted  values  a  useful  tool  in 
predicting  the  energy  usage  for  roof-spray  cooled  buildings. 

The  uncertainty  analysis  described  in  Appendix  B  indicate  that  the 
thermal  characteristics  of  the  composite  roof  are  important  factors  in  dry  and 
damp  roof  cooling  load  calculations  and  must  be  accurately  described  to 
ensure  the  accuracy  of  the  predicted  heat  gains  through  the  roof  .  Two  other 
variables  that  have  a  large  affect  on  the  uncertainty  of  the  predicted  heat  gain 
are  the  external  convective  heat  transfer  coefficient,  h„,  and  the  solar  radiation 
transferred  to  the  roof  surface  as  a  result  of  the  incident  solar  radiation  and 
the  absorptance  of  the  surface  of  the  roof  Since  h„  is  largely  a  function  of  the 
wind  speed  across  the  roof,  accurate  measurement  of  the  wind  speed  at  the 
building  location  is  important  to  allow  the  accurate  prediction  of  h,.  The 
effect  of  changes  in  the  absorptance  of  the  roof  surface  indicates  the 
importance  of  this  measurement  to  allow  the  accurate  prediction  of  solar 
radiation  transferred  to  the  roof  It  also  indicates  the  value  of  low 
absorptance  roof  surfaces  for  dry  and  roof-spray  cooled  buildings.  The  need 
to  accurately  measure  the  wind  speed  and  incident  solar  radiation  is  important, 
not  only  for  experimental  activities,  but  for  use  in  the  prediction  of  the 
required  application  of  water  to  the  roof  surface.  Therefore  the  installation  of 
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anemometric  and  heliopyrometric  devices  as  part  of  the  roof-spray  cooling 
system  may  be  desirable  for  the  optimization  of  roof  sprays  for  larger 
buildings. 


Recommendations 


The  following  is  a  list  of  suggested  improvements  that  can  be  included 
in  the  modeling  of  roof-spray  cooled  structures. 

1)  The  variation  in  wind  speed  over  the  roof  can  be  incorporated  into 
the  convective  heat  transfer  portion  of  the  outside  load  calculation. 

2)  The  effect  of  variations  in  cloud  cover  and  sky  clearness  on  the  net 
radiative  heat  exchange  between  the  roof  and  the  atmosphere  can  be  included 
in  the  outside  load  calculation. 

3)  The  predicted  evaporative  heat  transfer  rates  can  be  extended  to 
predict  the  variation  in  the  required  rate  of  delivery  of  water  to  the  roof  to 
optimize  the  cooling  effect  without  the  application  of  excess  water. 

4)  Various  roof  surface  materials  such  as  aggregate  ballast  of  different 
sizes  and  colors  can  be  compared  to  smooth  roof  surfaces  to  determine  any 
effects  on  the  heat  transfer  at  the  roof  surface. 

The  mathematical  formulation  provided  acceptable  results  for  buildings 
maintained  at  constant  inside  temperatures.  The  validation  of  this  model  for 
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buildings  that  are  maintained  at  different  temperatures  during  a  daily  or 
weekly  cycle  is  needed  because  of  the  common  usage  of  temperature  setbacks 
as  an  energy  conservation  measure. 

Any  further  experiments  on  roof-spray  cooling  should  include  a  study  of 
the  roofs  thermal  properties  including  surface  absorptivities  and  emissivities 
in  addition  to  ambient  weather  conditions  that  include  dry-bulb  temperature, 
humidity,  wind  speed,  pressure,  solar  radiation,  cloud  cover,  and  precipitation. 


APPENDIX  A 

COMPUTER  CODE  LISTING 


This  section  includes  the  BASIC  computer  code  that  was  used  to 
calculate  the  predicted  heat  fluxes  at  the  exterior  and  interior  surfaces  of  the 
roof  with  roof-spray  cooling  in  Bryan,  Texas.  Although  the  code  includes  the 
data  for  many  input  parameters  in  the  text,  the  code  could  be  easily  modified 
to  read  these  parameters  from  an  input  data  file  or  to  prompt  the  user  for  the 
input  values.  The  code  could  also  be  easily  translated  into  FORTRAN  if 
desired  for  compatibility  with  other  software  or  data  files. 

The  following  code  calculates  the  predicted  heat  fluxes  at  the  exterior 
and  interior  roof  surfaces  of  a  test  building  that  Somasundaram  et  al.  [3]  used 
for  roof-spray  cooling  experiments.  In  this  example  combined  convective  and 
radiative  heat  transfer  coefficients  were  used,  but  the  code  allows  the  use  of 
convective-only  heat  transfer  coefficients  if  the  radiative  heat  fluxes  are  not 
set  equal  to  zero.  The  external  convective  heat  transfer  coefficient,  h„,  can 
also  be  varied  if  needed  to  account  for  large  variations  in  wind  speed,  v. 

The  iterative  techniques  used  in  this  code  proved  to  be  effective  in  the 
required  simulation.  However,  the  speed  of  the  solution  could  be  improved 
through  the  application  of  more  sophisticated  bracketing  and  convergence 
techniques. 
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REM  ROOF-SPRAY  COOLED  ROOF  AT  BRYAN,  TEXAS  (IN  FILE 
ATMYBRTX.BAS) 

REM  HEAT  FLOW  CALCULATION  BY  TRANSFER  FUNCTION 
METHOD 

OPEN  "C  ATMBRTX  PRN"  FOR  OUTPUT  AS  #1 
WRITE  #1,  "Hour",  "OMEGA",  "Tset",  "Ta",  "Twet",  "TrO",  "Try", 
"Qsolar",  "Qevap",  "QcondO",  "Qe",  "QoutEXP",  "QinEXP", 
"qCONVo",  "pWR",  "pW",  "  " 

LET  T=0 

5  LET  H=T 

6  IF  H>=24  THEN  H=H-24 

IF  H>=24  GOTO  6 

DIM  Tr(-6  TO  120),  TE(-6  TO  120),  QE(-6  TO  120),  Try(-6  TO  120), 
1(0  TO  120),  QEXPout(0  TO  120),  QEXPin(0  TO  120) 

REM  Tdry-bulb,  Twet-bulb,  PRESSURE,  HUMIDITY,  AND  SOLAR 
VALUES  FROM  TYPICAL  METEOROLOGICAL  YEAR 
WEATHER  DATA  AVERAGED  FOR  PERIOD  17-30  JULY  FOR 
WACO,  TEXAS 

REM  WEATHER  FILES  OBTAINED  FROM  ACROSOFT 

INTERNATIONAL,  INC.,  DENVER,  CO  AND  PROCESSED  BY 
MICRO-DOE2  SOFTWARE 

REM  EXPERIMENTAL  DATA  FROM  SOMASUNDARAM  et  al.  [1988]  pp 
1103-1104 

IF  H=0  THEN  TW=69:  TAMB=79.3:  P=29.45;  W=  031;  QSOL=0; 
QEXPout=-3;  QEXPin=.45 

IF  H=1  THEN  TW=69;  TAMB=77.8:  P=29.47;  W=  0136;  QSOL=0; 
QEXPout=-3;  QEXPm=.34 

IF  H=2  THEN  TW=68.9;  TAMB=76.5;  P=29.47;  W=.0138;  QSOL=0; 
QEXPout=-3.15;  QEXPin=.27 

IF  H=3  THEN  TW=68.7;  TAMB=75.3;  P=29.47;  W=.014;  QSOL=0; 
QEXPout=-3.15;  QEXPin=.24 

IF  H=4  THEN  TW=68.6.  TAMB=74.6;  P=29.48;  W=.014;  QSOL=0; 
QEXPout=-3.15;  QEXPin=  21 

IF  H=5  THEN  TW=68.3;  TAMB=73.7;  P=29.48;  W=.014;  QSOL=0; 
QEXPout=-3.15;  QEXPin=,17 

IF  H=6  THEN  TW=68.1;  TAMB=72.9;  P=29.49;  W=  014;  QSOL=0; 
QEXPout=-3.15;  QEXPin=  14 

IF  H=7  THEN  TW=69.7;  TAMB=76;  P=29.51;  W=.0145;  QSOL=15.5; 

QEXPout=-2.47;  QEXPin=-  21 
IF  H-8  THEN  TW=71.1;  TAMB=79.1;  P=29.51;  W=.0149; 
QSOL=104.9;  QEXPout=- 1 . 1 ;  QEXPin=-  2 
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IF  H=9  THEN  TW=72  6:  TAMB=82  2:  P=29  51  W=  0153 
QSOL=170  3:  QEXPout=l  64:  QEXPin=-  24 
IF  H=10  THEN  TW=72.7;  TAMB=84  9:  P=29  51  W=  0149 
QSOL=202.6  QEXPout=3  83:  QEXPin=-  24 
IF  H=ll  THEN  TW=73  1  TAMB=87  4:  P=29.51:  W=  0146: 

QSOL=215.4:  QEXPout=5  34:  QEXPin=-  07 
IF  H=12  THEN  TW=73  3;  TAMB=90  1:  P=29  49:  W=  0141 
QSOL=218  2:  QEXPout=6  59:  QEXPin=  09 
IF  H=13  THEN  TW=73  1:  TAMB=91  3:  P=29  49:  W=  0136: 

QSOL=213.1:  QEXPout=7  12:  QEXPin=  38 
IF  H=14  THEN  TW=72.9;  TAMB=92.7:  P=29  47:  W=  0132: 

QSOL=206.9:  QEXPout=5.62:  QEXPin=  48 
IF  H=15  THEN  TW=72.9;  TAMB=94:  P=29  43:  W=  0129: 

QSOL=205.4:  QEXPout=5.2:  QEXPin=  75 
IF  H=16  THEN  TW=72.5;  TAMB=93.3;  P=29.41:  W=  0128: 

QSOL=  199.1;  QEXPout=3.7:  QEXPin=l  09 
IF  H=17  THEN  TW=72.1;  TAMB=92.7:  P=29.41;  W=  0126: 

QSOL=  178.7:  QEXPout=3.3;  QEXPin=l  06 
IF  H=18  THEN  TW=71.5;  TAMB=92:  P=29.41;  W=  0122: 

QSOL=124.2:  QEXPout=1.64.  QEXPin=1.06 
IF  H=19  THEN  TW=70.8.  TAMB=88.7:  P=29.41:  W=.0124: 

QSOL=4.7:  QEXPout=-.13:  QEXPin=1.35 
IF  H=20  THEN  TW=70.3:  TAMB=85.9:  P=29.42:  W=.0127;  QSOL=0: 
QEXPout=- 1.1:  QEXPin=  1.16 

IF  H=21  THEN  TW=69.8:  TAMB=82.8:  P=29.43:  W=.0129:  QSOL=0: 
QEXPout=-2.33:  QEXPin=.99 

IF  H=22  THEN  TW=69.5:  TAMB=81.7:  P=29.44:  W=.0131:  QSOL=0: 
QEXPout=-2.74:  QEXPin=.75 

IF  H=23  THEN  TW=69.1:  TAMB=80.4:  P=29.45:  W=.0130:  QSOL=0: 
QEXPout=-2.86:  QEXPin=.49 

IF  H=24  THEN  TW=69:  TAMB=79.3:  P=29.45:  W=.031:  QSOL=0: 

QEXPout=-3:  QEXPin=.45 
REM  END  OF  TMY  DATA 

P=P*3386  4:  REM  CONVERT  PRESSURE  FROM  INCHES  HG  TO  Pa 
Pw=(P*W)/(.62198+W):  REM  CALCULATE  PARTIAL  PRESSURE  OF 
WATER  IN  MOIST  AMBIENT  AIR 
LET  QSOL=.6*QSOL*3  155:  REM  APPLY  ROOF  ABSORPTANCE 
AND  CONVERT  FROM  Btuh/sf  TO  W/m**2 
Twet=(TW-32)/l. 8+273. 15;  REM  CONVERT  WET  BULB  TEMP 
FROM  F  TO  K 

TAMB=(TAMB-32)/1.8+273.15:  REM  CONVERT  DRY  BULB  TEMP 
FROM  F  TO  K 
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REM  FROM  ROOF  GROUP  #166  BUILD  ON  TRANE  TRACE  600 

ROOF  GENERATOR  UTILITY  BASED  ON  SOMASUNDARAM 
et  al  TEST  BUILDING 

REM  DESCRIPTION  OF  ROOF,  CONDUCTION  TRANSFER 
FUNCTION  COEFFICIENTS 

BO  =  4  40399E-09*5  6783;  REM  CONVERT  Bn  AND  Cn  TO  SI  UNIT 
B1  =  1.461269E-05*5.6783;  B2  =  2.489753E-04*5.6783 
B3  =  5.109677E-04*5  6783;  B4  =  2.08079E-04*5.6783 
B5  =  1.789627E-05*5.6783;  B6  =  2.828292E-07*5.6783 
Cn  =  001*5.6783:  DO  =  1 

DI  =  -1.974663:  D2  =  1  273348:  D3  =  -.3056189 
D4  =  .02383278;  D5  =  -4.312459E-04:  D6  =  1.107164E-06 
REM  TIME  LAG  =  9  HOURS,  Utable  =  .061*5.6783,  W/(m**2*K) 
REM  DF  =  227,  DECREMENT  FACTOR,  AMPLITUDE  REDUCTION 
FACTOR  FROM  OUTSIDE  HEAT  FLUX  TO  INTERIOR  LOAD 
Ho  =  22  7:  REM  EXT  CONVECTIVE  AND  RADIATIVE  HEAT 
TRANSFER  COEFFICIENT,  W/(m**2*K) 

Hi  =  8  35:  REM  INT  CONVECTIVE  AND  RADIATIVE  HEAT 
TRANSFER  COEFFICIENT,  W/(m**2*K) 

EPS  =  85:  REM  EMISSIVITY  OF  ROOF 
SIG  =  5  67E-08.  REM  STEFAN-BOLTZMAN  CONSTANT, 
W/(m**2*K**4) 

Rr  =  2  079:  REM  THERMAL  RESISTANCE  OF  ROOF,  m**2*K/W 
TSET  =  299  81:  REM  AMBIENT  TEMPERATURE  OF  INSIDE  AIR 
PLENUM,  K 

LET  Try(-l)  =  280:  REM  PROVIDES  STARTING  TRIAL  TEMP  FOR 
HOURO 

IF  T>=24  THEN  Try(T)=Try(T-24)-  1.0001:  GOTO  20 
Try(T)  =  Try(T-l)  -  20 

20  Fry  =  Try(T)*  18:  CONVERT  TEMP  FROM  K  TO  F 

REM  CALCULATE  SATURATION  PRESSURE  OF  WATER  ON 
ROOF  SURFACE  FROM  MATHUR  [1988] 

PARTO  =  670.2012:  PARTI  =  -5.325521  *Fry 

PART2  =  .0159464*Fry**2:  PART3  =  -.00002134061  *Fry**3 

PART4  =  1.077853E-08*Fry**4 

Pwr  =  (PART0+PART1+PART2+PART3+PART4)*6894.8:  REM  Pa 
QEVAP  =  .073814*(Pwr-Pw):  REM  CALCULATE  EVAPORATIVE 
HEAT  TRANSFER  DUE  TO  PRESSURE  DIFFERENCE  OF 
WATER  VAPOR  AT  ROOF  SURFACE  AND  IN  AMBIENT  AIR 
REM  Tr(T)  =  TAMB+Rr*(-QEVAP-QRAD-QCONV+QSOL) 

QRADo  =  EPS*SIG*(Try(T)**4-(TAMB-3)**4) 

QCONVo  =  Ho*(Try(T)-TAMB) 


37 


Tri  =  Try(T)+Rr*(QRADo+QEVAP+QCONVo-QSOL)  REM  CALC 
INT  ROOF  TEMP,  K 

QCONDo  =  -QRADO-QEVAP-QCONVo+QSOL 
QCONVi  =  Hi*(Tri-TSET) 

QRADi  =  EPS*SIG*(Tri**4-TSET**4) 

REM  QCONDi  =  QCONVi=QRADi:  REM  INT  SURFACE  ENERGY 
BALANCE 

Tr(T)  =  Tri+Rr*(QCONVi+QRADi) 

IF  ABS(Tr(T)-Try(T))<  1  GOTO  40:  REM  CHECK  FOR 
CONVERGENCE  OF  TRIAL  ROOF  TEMP  WITH 
CALCULATED  TEMPERATURE 

REM  SET  STEP  SIZE  TO  CALCULATE  OUTSIDE  ROOF  TEMP 
IF  Tr(T)>285  AND  Tri>0  THEN  I(T)=  0001:  GOTO  30 
IF  Tr(T)>270  AND  Tri>0  THEN  I(T)=  001:  GOTO  30 
IF  Tr(T)>0  AND  Tri>0  THEN  I(T)=  01:  GOTO  30 
IF  Tr(T)<0  OR  Tri<0  THEN  I(T)=  1 

30  LET  Try(T)  =  Try(T)+I(T):  REM  ADD  ITERATIVE  STEP  TO 
PREVIOUS  TRIAL  ROOF  TEMPERATURE 
IF  Try(T)>Tr(T-l)+20  THEN  PRINT  "ERROR  ON  ITERATIVE 
LOOP":  GOTO  50 
GOTO  20 

40  REM  SOL-AIR  TEMPERATURE  WILL  BE  TAKEN  AS 

CALCULATED  ROOF  SURFACE  TEMPERATURE 
LET  TE(T)  =  Tr(T)-273.I5:  REM  CONVERT  ROOF  TEMP  FROM  K 
TO  C 

REM  CALCULATE  PREDICTED  LOAD  ON  INTERIOR  FROM  ROOF 
BY  TRANSFER  FUNCTION  METHOD 
QEI=B0*TE(T)+Bl*TE(T-l)+B2*TE(T-2)+B3*TE(T-3) 
QE3=B4*TE(T-4)+B5*TE(T-5)+B6*TE(T-6) 

QE3=D  1  *QE(T- 1  )+D2*QE(T-2)+D3  •QE(T-3) 
QE4=D4*QE(T-4)+D5*QE(T-5)+D6*QE(T-6)+(TSET-273.15)*Cn 
QE(T)  =  QE1+QE2-QE3-QE4 
REM  WRITE  OUTPUT  DATA  TO  FILE 

WRITE  #1,  T,  W,  TSET,  TAMB,  TWET,  Tr(T),  QSOL,  QEVAP, 

QCONDo,  QE(T),  QEXPout(T),  QEXPin(T),  QCONVo,  Pwr,  Pw 
LET  T  =  T+1:  REM  START  NEXT  HOUR  SIMULATION 
IF  T>120.5  GOTO  50:  REM  STOP  SIMULATION  AFTER  120  HOURS 
GOTO  5 

50  REM  END  OF  WHILE  LOOP 

CLOSE  #1:  REM  CLOSE  OUTPUT  FILE  ATMYBRTX.PRN 
END 
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REM  *********nOMENCLATURE********** ******************* 

Bn  =  CONDUCTION  TRANSFER  FUNCTION  COEFFICIENTS 

Cn  =  CONDUCTION  TRANSFER  FUNCTION  COEFFICIENTS 

Dn  =  CONDUCTION  TRANSFER  FUNCTION  COEFFICIENTS 

EPS  =  EMISSIVITY 

Fry  =  TRIAL  ROOF  TEMPERATURE,  DEGREES  R 

H  =  HOUR  OF  DAY,  MODULUS  24  CLOCK,  HOURS 

Hi  =  INTERNAL  CONVECTIVE(/  RADIATIVE)  HEAT  TRANSFER 

COEFFICIENT,  W/(M**2*K) 

Ho  =  EXTERNAL  CONVECTIVE(/  RADIATIVE)  HEAT  TRANSFER 
COEFFICIENT,  W/(M**2*K) 

P  =  TOTAL  PRESSURE  OF  MOIST  AIR,  INCHES  HG,  PSIA  OR  Pa 

Pw  =  PARTIAL  PRESSURE  OF  WATER  IN  MOIST  AIR,  Pa 

Pwr  =  PARTIAL  PRESSURE  OF  WATER  AT  ROOF  SURFACE  (Psat 

AT  Troof),  Pa 

Rr  =  THERMAL  RESISTANCE  OF  ROOF,  m**2*K/W 

QCONDo  =  CONDUCTIVE  HEAT  FLUX  INTO  ROOF  FROM  OUTER 
SURFACE,  W/m**2 

QCONVi  =  CONVECTIVE  HEAT  FLUX  AT  BOTTOM  OF  ROOF  TO 
ROOM,  W/m**2 

QCONVo  =  CONVECTIVE  HEAT  FLUX  AT  TOP  ROOF  TO  AMBIENT, 
W/m**2 

QE(T)  =  PREDICTED  HEAT  FLUX  INTO  ROOM,  W/in**2 
QEVAP  =  EVAPORATIVE  HEAT  FLUX  FROM  MOIST  ROOF,  W/m**2 
QRADi  =  RADIATIVE  HEAT  TRANSFER  FROM  ROOF  INNER 
SURFACE  TO  AMBIENT  ROOM,  W/m**2 
QRADo  =  RADIATIVE  HEAT  TRANSFER  FROM  ROOF  OUTER 
SURFACE  TO  AMBIENT,  W/m**2 

QSOL  =  SOLAR  RADIATION,  INCIDENT  AND  ABSORBED,  W/m**2 
SIG  =  SIGMA-STEFAN-BOLTZMAN  CONSTANT,  5  6697E-08 
W/(ra**2*K**4) 

T  =  NUMBER  OF  HOURS  STARTING  AT  MIDNIGHT  ON  FIRST 

SIMULATION  DAY,  H 

TAMB  =  AMBIENT  DRY  BULB  TEMPERATURE,  F,  C  OR  K 
TE(T)  =  CALCULATED  EQUIVALENT  ROOF  TEMPERATURE  =(Tr),  C 
Tr(T)  =  ROOF  SURFACE  TEMPERATURE  AT  HOUR  T,  C  OR  K 
Tri  =  CALCULATED  INSIDE  ROOF  SURFACE  TEMPERATURE,  K 
Try(T)  =  TRIAL  ROOF  SURFACE  TEMPERATURE,  K 
TSET  =  AMBIENT  ROOM/ AIR  PLENUM  TEMPERATURE,  C  OR  K 

TW  =  AMBIENT  WET  BULB  TEMPERATURE,  F,  C  OR  K 
W  =  OMEGA-HUMIDITY  RATIO,  MASS  OF  WATER  PER  UNIT 

MASS  OF  DRY  AIR 


APPENDIX  B 

UNCERTAINTY  ANALYSIS 


To  estimate  the  accuracy  of  the  predicted  results  it  is  necessary  to 
quantify  the  uncertainty  of  the  individual  variables  and  its  effect  on  the 
uncertainty  of  the  model  results.  In  this  thesis  the  result  is  the  prediction  of 
the  heat  gain  through  the  roof,  q,.  The  variables  include  ambient  temperature, 
humidity  ratio,  solar  heat  flux  at  the  roof  surface,  inside  and  outside 
convective  heat  transfer  coefficients,  thermal  resistance  of  the  roof,  and  net 
radiative  heat  exchange  with  the  sky. 

The  uncertainty  analysis  method  chosen  is  the  sequential  perturbation  of 
the  data  reduction  program  described  in  Appendix  A.  This  method  is 
described  in  detail  by  Moffat  [19],  This  method  is  based  on  determining  the 
effect  of  each  variable's  uncertainty  on  the  baseline  predicted  result.  After  the 
result  of  the  baseline  data  set  are  calculated  results  are  determined  "once  more 
for  each  variable,  with  the  value  of  the  variable  increased  by  its  uncertainty 
interval  (and  all  other  [variables]  returned  to  their  baseline  values)".  [19]  The 
difference  between  the  returned  perturbed  results  and  the  baseline  represents 
the  contribution  of  each  variable’s  uncertainty  on  the  uncertainty  of  the  result's 
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overall  uncertainty  The  uncertainty  of  the  result  is  determined  by  squaring 
each  contribution,  summing,  and  taking  the  square  root. 

The  expression  used  to  calculate  the  uncertainty  W  of  a  function  is 


i 

W={|(Wx.)Y  (B-1) 


where  is  any  of  the  variables  which  quantify  the  function. 

The  uncertainty  interval  for  each  variable  has  been  estimated  based  on 
Somasundaram  et  al.  [5]  and  Houghten  et  al.[I]  .  The  results  of  the  sequential 
perturbation  for  the  roof-spray  cooled  concrete  roof  at  Pittsburgh  are 
provided  in  Table  B.  1 . 


Table  B.l  Uncertainty  Analysis  by  Sequential  Perturbation  for  Concrete 
Roof  at  Pittsburgh 


Case/ 

Uncertainty 

Perturbed 

Difference  From 

Perturbed 

Int^al 

Result 

Baseline 

Variable 

(W/m") 

(W/m^) 

Basdine 

-4.73 

T. 

0.5  K 

-4.28 

0.45 

u 

0.001  kg,^g*y^ 

-4.72 

0.01 

Q-.-  and  a 

lOW/m^ 

-4.27 

0.46 

h, 

2  W/(m^  K) 

-4.71 

0.02 

K 

4  W/(m'  K) 

-4.28 

0.45 

R, 

0.05  m'  K/W 

-4.42 

0.31 

T 

^•et 

0.5  K 

-5.42 

0.69 

SR 

10  W/m' 

-5.13 

0.39 
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The  calculation  of  the  total  uncertainty,  W,  is  shown  by 
W=  {(0,45)'+(0.01)'+(0.46)'+(0.02)'+(0  45)'+(0.31)'+(0.69)'+(0  39)'}“’ 
which  yields  an  uncertainty,  W,  of  12  W/m'  Similar  calculations  for  the 
other  roofs  yielded  uncertainties  of  0.2  W/m'  for  the  predicted  0.40  W/m' 
predicted  average  roof  heat  gain  for  the  spray  cooled  roof  at  Bryan,  Texas,  12 
W/m'  for  the  predicted  2.4  W/m'  predicted  average  for  the  dry  roof  at  Bryan, 
and  1.4  W/m'  for  the  predicted  -4.2  W/m'  average  for  the  spray  cooled  pine 
plank  roof  at  Pittsburgh. 
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